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Abstract 



Context. A planetary nebula (PN) candidate was discovered during FORS imaging of the Local Group dwarf galaxy Phoenix. 
Aims. Use this PN to complement abundances from red-giant stars. 

Methods. FORS spectroscopy was used to confirm the PN classification. Empirical methods and photoionization modelling were 
used to derive elemental abundances from the emission line fluxes and to characterize the central star 

Results. For the elements deemed most reliable for measuring the metallicity of the interstellar medium (ISM) from which the PN 

formed, [O/H] 0.46 and [Ar/H] 1.03. [0/H] has lower measurement errors but greater uncertainties due to the um'esolved issue 

of oxygen enrichment in the PN precursor star. 

Conclusions. Earlier than 2 Gyr ago (the lower limit of the derived age for the central star) the ISM had Z = 0.002-0.008, a range 
slightly more metal-rich than the one provided by stars. Comparing our PN-to-stellar values to surveys for other dwarf Local Group 
galaxies, Phoenix appears as an outlier 

Key words, galaxies: dwarf - galaxies: individual: Phoenix - planetary nebulae: individual: PN Phoenix 101:51:05.46-44:26:55.28 



1. Introduction 

1.1. The Phoenix dwarf and its metallicity 

Phoenix is a Local Group dwarf galaxy located 450 kpc from the 
Milky Way [My = -10.1, Ly = 9 x 10^ Lq, M= 3.3 x IO^'Mq, 
(m - M)o - 23.2; Mateo 11998] ). It has little curren t star for- 
mation (Martmez-Delgado, Gallart & Aparicio |1999l hereafter 
M99; Held, Saviane & Momany T999, hereafter H99) a nd lit- 
tle or no H I gas (St-Ge rmain et al. il999i G allart et al. 120011 
Irwin & Tolstoy I2002t Young et al. 12007b . which warrants 
its classification as a transition-type dwarf galaxy. From the 
color and width of the red-giant branch (RGB) one can infer 
quite a low metallicity for its intermediate-age and old stars: 
[Fe/H]= -1.85, o-([Fe/H])^ 0.5 dex (H99). However, metal- 
licity estimates from the position of the RGB alone are quite 
uncertain in the case of composite stellar populations due to 
the resulting age-metallicity degeneracy (see, e.g. Pont et al. 
I2004I I. With a color-magnitude diagram (CMD) reaching the 
oldest main-sequence turnoffs, the RGB age-metallicity degen- 
eracy may be partially broken since ages can be derived from 
the main-sequence, and Z(t) can be estimated from the global fit 
of the distribution of stars in the CMD using stellar evolution 
models. In this way. Hidalgo, Aparicio & Martmez-Delgado- 
Delgado (120071 H07) have found that the current metallicity 
should be about Z = 0.0015. 
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In any case, spectroscopic data provide a more direct mea- 
surement of the metallicity and metallicity distribution of the 
stars in a galaxy. In the case of Phoenix, Ca n triplet spectra, 
which provide global metallicities, could be obtained for the 
brightest RGB stars in the galaxy. Measuring abundances di- 
rectly, through high-resolution spectroscopy of individual stars, 
is not feasible due to the large distance of the galaxy. Moreover, 
there are no known Hii regions, so the only other opportu- 
nity to obtain information on the abundances of particular el- 
ements is through a PN. During a FORSl campaign (63.1- 
0642), a PN candidate was identified near the center of the 
galaxy, and in this paper we present our follow-up spectroscopy. 
The object was discovered due to its relative brightness in Ha 
as compared with the V band, together with its star-like ap- 
pearance (see Fig. [Til. Following current convention we name 
this object PN Phoenix JOl :5 1 :05 .45-44:26:55.28 but, in a light- 
hearted spirit, we choose to call it Bennu (the ancient Egyptian 
name for the phoenix). 



1.2. Chemical abundances from PNe 

PNe are the post-AGB (asymptotic giant branch) stage of low- 
to intermediate-mass stars. Their chemical abundances are those 
of the envelope of the precursor AGB star which was ejected 
via a stellar wind. The composition of the nebulae are a mix- 
ture of elements affected by the preceding nucleosynthesis and 
dredge-up cycles (e.g. He, C, and N) and elements unadulterated 
by stellar evolution (e.g. Ar, S, and in most cases also O and 
Ne). Measuring PN abundances is useful since they complement 
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Figure 1. Finding chart for the Phoenix planetary nebula, ob- 
tained from FORSl V-band imaging. The lower two insets show 
the appearance of the nebula in the Ha band (left) and V band 
(right). North is up and east is to the left. 



abundance determinations from stars, and contribute to the study 
of galactic chemical evolution (e.g.. Shields |2002| |. For Phoenix, 
obtaining PN abundances is especially important, as this is one 
of the most metal poor Local Group galaxies. While with only 
one PN we are dealing with low number statistics, we note that 
the scatter about the mean oxygen abundance for PNe (the most 
often considered element in PN-galactic evolution studies) in 
other galaxies is usually less than a factor of two (Stasihska & 
Izotov 120031). Lacking any other determination, a sole PN can 
still provide crucial information. 



2. Observations and data reduction 

The spectroscopic observations were carried out in service mode 
on the night of November 26, 2005. The target can be identified 
in Fig. [H which shows part of one FORS2 frame taken in the 
V band. Three FORSl spectra of 1730s each were taken, with 
the spectrograph in long-slit mode, a CK'? slit, and with grism 
300V+10. No order sorting filter was used, in order to have a 
higher throughput in the blue part of the spectrum, and the re- 
sponse function was then determined by observing a blue stan- 
dard star in the same configuration. Some second-order contam- 
ination is expected at A > 6600 A, but it is less than 4% for 
the [S ii] lines at 6725 A. The nebula was observed at an air- 
mass ~ L04 and under a seeing of ~ 1", while approaching and 
crossing the meridian. During the evening twilight, the DA spec- 
tropho tome trie standard EG21 from Hamuy et al. (I1994I I was 
observed. The instrument mode was slightly different than the 
one used for the science target: a 5" slit was created by placing 
the 19 movable slits of the MOS mode side by side, and align- 
ing them along the parallactic angle; moreover, only the central 
part of the CCD was read out. The science data were instead 



taken with one of the fixed-width slit masks, and the full CCD 
was read out. The two sets of calibration frames (biases, arc, and 
flat-fields) for the two configurations were taken in the morning 
following the observations. 

Reductions were carried out using a customized version of 
the EFOSC2 quick-look toofl, which is based on the ESO- 
MIDAS data reduction system. The usual bias subtraction and 
flat-fielding (from quartz lamps) were performed. He-i-HgCd+Ar 
arc frames were used to compute the 2-D wavelength calibration. 
The spectra were then linearly rebinned with a constant step of 
A/l = 2.64 A, with a resulting wavelength coverage of 3707- 
8289 A at a resolution of ~ 8.2 A FWHM. 

The two-dimensional sky frame was created by first sam- 
pling the sky spectrum in two windows flanking the target spec- 
trum, and then fitting the spatial gradient with a one-degree 
polynomial. After sky-subtracting, the extraction window for the 
emission-line spectrum was chosen based on the spatial profile 
of the Ha line, adjusting in order to reach an optimum signal- 
to-noise ratio (SNR). The flux of the spectrophotometric stan- 
dard star was instead summed over almost the entire spatial pro- 
file, leaving out just the two sky windows. To compute the re- 
sponse function, the instrumental spectrum was corrected for at- 
mospheric extinction, divided by the exposure time, and then 
divided by the published spectrum, and the ratio was fit with a 
12* de gree polynomial. Note that the different instrument setup 
of the standard does not affect the response function, which 
can therefore be applied to the science spectra. The three pre- 
processed, wavelength-calibrated, sky-subtracted and extracted 
spectra of the PN were flux-calibrated using this response func- 
tion. Finally, the fiducial PN spectrum was obtained as a median 
of the three individual exposures, and it is shown in Fig.|2] The 
typical SNR measured on the continuum is five or better, while 
the [O III] /14363 Une was detected with a SNR of ~ 24. The con- 
tinuum is due in part to the nebula and in part to the unresolved 
light of the Phoenix stellar population. Photoionization models 
(see below) predict that 39% of the observed continuum inten- 
sity at 5000 A is due to the nebula. Because of the low counts, 
no absorption lines can be detected, and no slope is visible. 



3. Spectral measurements 

Emission line fluxes and centroids were measured with the ELF 
Gaussian fitting routines of DIPSO (Howarth et al. 1 19981 1 and 
independently checked with the fitting routines of the MIDAS 
package. The FWHM for faint lines were fixed from near-by 
bright lines and were a free parameter for unresolved blends. The 
fluxes are reported in Table[T| The measured blended H/3+He ii 
flux is 1.98x10 '^ergcm^^s' and the adopted de-reddened and 
de-blended Hy6 intensit}0 is 3.50xlO""'ergcm"^ s"'. This value, 
which is used in the models of Section 4.2, incorporates a cor- 
rection of +28% as the employed slit had a width of 0."7 and the 
PN was observed under ~l" seeing. 

The continuum in the very blue and red is fairly noisy, and 
a conservative approach was taken to deciding between a true 
line and a noise spike. This is especially important for the Hei 
lines - because of the dependence of the calculated oxygen abun- 
dance on the helium abundances when using the TCF' method 
(see below) - and the [S ii] ratio, which is used to calculate the 
electron density. The sky is noisy in the [S ii] region, and the 
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' Correction for the He ii contribution is calculated from theoretical 
He II line ratios to the uncontaminated He ii 4686 A flux. 
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Figure 2. The calibrated and sky-subtracted spectrum of Bennu, above, and again below, scaled to highlight the fainter lines. 



[S n] fluxes very low. To estimate the effect this has on the mea- 
surement of the [S n] ratio we performed a sky subtraction in the 
standard manner (using the sky spectrum extracted from the slit 
and resulting in the fluxes quoted in Table 1) and a sky subtrac- 
tion using rather a polynomial fit to the sky in this spectral re- 
gion. The [S ii] ratio from the poly-fit result is up to 70 % greater 
that from the true-sky subtraction. Regarding He+, only the He i 
/15876 line, de-blended from He ii /15869.02 (5-29), is consid- 
ered in our analysis as the 4471 A and 6678 A lines are upper 
limit detections. 



The Ha/H/S Balmer line ratio was used to calculate the ex- 
tinction constant, yielding a value of c(Hy8)=0.16. The Hen 
Pickering series is a contaminant of the H i lines in high exci- 
tation nebulae such as this one and their contribution was esti- 
mated via their theoretical ratios to He ii i4686 following Storey 
& Hummer ( fT995] l; 13.5% of the /I4686 flux at Ho- is due to He ii 
Pi (4-6) and 5.1 per cent at H/S is due to Pi (4-8), under case 
B conditions. The intrinsic value of the Hi ratio to compare to 
when computing c(Hy6) depends on the electron temperature, T^., 
and density, n^; we iterated once between determining Tg and «e 
before computing c(Hy8). As the Hen /1/14541.59, 5411.53 lines 
are well detected, we also computed, to compare, the redden- 
ing using the Hen /i4686//l5411 and /15411//I4542 ratios; this 
yielded a weighted average of 0.19 + 50%. We adopted the Hi 
Balmer value to deredden the observed line fluxes, and used the 
Galactic extinction law of Howarth (1983). 



A 16% correction has been applied to the [O ii] A3127 dou- 
blet which is blended with the H 14 ^3721.9, [S m] /13721.6 and 
H 13 /13734.4 lines: the [S iii] /13721 flux was estimated from a 
comparison with the de-reddened flux of [S iii] /16312.1, which 
originates from the same upper level. The latter line is blended 
with He II /i6310.8 (5-16) in high excitation PNe and that flux 
was retrieved via its theoretical ratio relative to He ii /14686. The 
/(H 13)//(H;e) and /(H 14)//(Hy6) intensity ratios were estimated 
using theoretical Hi line emissivities as above. 

In Table [U suggested ion identifications are given. The bluest 
Hi line ratios, at A < 3968, are a factor of 1.3-1.8 in excess 
of that expected theoretically, suggesting that the adopted blue 
continuum, possibly aff'ected by interstellar absorption lines, is 
too low. H I /13771 is not detected in the spectrum. 



4. Chemical abundances and central star 
parameters 

4. 1 . Empirical analysis 

A determination of chemical abundances was first carried out us- 
ing the semi-empirical ICF (ionization correction factor) scheme 
reported in Kingsburgh & Barlow (1994). The abundances of 
each observed ionic species relative to H^ are added up and 
multiplied by their corresponding ICF to yield elemental abun- 
dances relative to hydrogen. Input values in these calculations 
are the plasma Te and rig, whose values were computed from the 



Table 1. Observed line fluxes. 







1 y/Lj 


ErT.(%) 


ID 


3727 


54.1 


26.83 


10 


[On] 


3749 


10.4 


5.73 


50 


Hi 


3798 


13.4 


7.37 


30 


Hi 


3836 


18.6 


10.20 


15 


Hi 


3869 


48.6 


28.05 


10 


[Ne in] 


3889 


24.5 


13.54 


10 


Hi, He I 


3969 


48.5 


26.93 


10 


Hi, [Nem] 


4069 


6.43 


3.66 


25 


[S 11]+^ 


4102 


51.4 


27.75 


10 


Hi 


4341 


95.1 


50.43 


10 


Hi 


4363 


15.0 


8.32 


10 


[O III] 


4541 


7.46 


4.07 


40 


Hen 


4657 


5.59 


3.02 


50 


C IV, [Fe III] 


4686 


182. 


98.06 


10 


Hen 


4714 


10.1 


5.44 


30 


[Ar iv]+'' 


4740 


5.06 


2.71 


35 


[Ar iv] 


4861 


198. 


100 


10 


Hi 


4959 


210. 


110 


5 


[O III] 


5007 


651. 


341 


5 


[Om] 


5412 


17.2 


8.68 


20 


Hen 


5876 


5.70 


2.77 


30 


He I 


6312 


3.11 


1.11 


50 


[S in] 


6563 


630. 


284 


5 


Hi 


6583 


19.9 


9.37 


15 


[Nil] 


6717 


3.57 


1.67 


50 


[S II] 


6731 


5.61 


2.62 


50 


[S II] 


7136 


5.95 


2.73 


25 


[Ar in] 



["] In units of 10"'^ erg cm"* s"'. 

[*] De-reddened and ratioed to IIy6 = 100. All H i lines de-blended from 
He n; [O n] 3727 de-blended from [S m] and two Hi lines. 

[^ Blends: [S n] 4069 + [Fe v] 4071 + C m 4069; [Ar iv] 471 1 + He i 
4713 + [Ne iv] 4714-^16. 



Table 2. Empirical ionic abundances" and ICF values for two 
electron densities. 



Ion 


500 


3700 


Ion 


500 


3700 


[On] 


17.1 


25.7 


[Sn] 


0.38 


0.54 


[O m] 


280. 


279. 


[S ra] 


3.87 


4.30 


icf(0) 


3.96 


3.96 


icf(S) 


1.82 


1.61 


[Nil] 


5.42 


5.62 


[Ar ra] 


0.93 


0.93 


icf(N) 


73.1 


50.0 


[Ar IV] 


1.82 


1.64 


[Ne ra] 


57.4 


57.4 


icf(Ar) 


1.01 


1.02 


icf(Ne) 


4.18 


4.30 









["] In 10 ^ by number relative to H. 



de-reddened [Om] /14363/(^5007 +^4959) and [S n] /167 17/6731 
ratios respectively. 

The same atomic data references as in Exter et al. ( I2004I I 
were used, with the exception of using helium where effective re- 
combination coefficients were taken from Smits ( 119961 ) for Case 
B and adopting corrections for collisional excitation contribu- 
tions as in Benjamin, Skillman and Smits (119991) . 

For the density measured from the [S ii] ratio: to show the 
effect that measurement uncertainties (in particular, sky subtrac- 
tion) has on the derived density we carried out here calcula- 
tions using the [S ii] ratio measured from the two sky-subtraction 
methods mentioned in Sec. 3. A higher density comes from the 



Table 3. Abundances of Bennu vs solar^ 



Element 


ne= 500 


= 3700 


Model 


Sun 


He/H 


.125(20) 


.125(20) 


.108(10) 


.098 


O/HxlO" 


1.18(14) 


1.21(17) 


1.7(13) 


4.9 


N/OxlO 


3.4(24) 


2.3(26) 


1.5(15) 


1.6 


Ne/OxlO 


2.0(30) 


2.0(26) 


1.4(15) 


1.8 


S/OxlO^ 


6.5(50) 


6.5(50) 


10.'' 


37. 


Ar/OxlO^ 


2.4(30) 


2.2(30) 


2.3(25) 


8.5 



[^] Proto-solar abundances from Asplund et al. ( I2005t . except for O/H 
from Allende-Prieto et al. l l2001t . Percentage errors in parenthesis. 

[*"] S/Ar solar; S/OxlO^^ may as well be 20 according to model out- 
puts. 



lower ratio (from the fluxes we adopted and list in Table 1) and 
a lower density from the higher ratio. The effect on the Te value 
is minor, 7^.= 16 600 + 1000 K, while changes from 500 to 
3700 cm"^. In Tables 2 and 3 we give the ionic and total abun- 
dances with this range of «e and values. The percentage 
Icr errors given in parenthesis in Table 3 include the flux and 
Tf. measurement errors and not those inherent to the empirical 
methods (see Exter et al. 120041 ). In Table 3, errors for the model 
abundances (see Sect. 4.2) correspond either to the scatter among 
several model fits or to the flux errors. 

Since the He n flux is very high, the nebula cannot be radi- 
ation bounded. Then the formulation of Kaler & Jacoby ( 119891) 
for estimaging the Teff of the central star only provides an upper 
limit to Teff (< 322 kK) and the one of Zijlstra & Pottasch {1^^ 
a lower limit to Lboi (> 270 L©). The method of Zanstra (119271 ) to 
obtain a reff(He ii) requires knowledge of the stellar continuum 
(see below), while the energy balance method of Stoy ( 119331 ) re- 
quires knowledge of important gas coolants outside the optical 
range. Self-consistent photoionization modelling based on the 
optical line spectrum is the only available method for the deter- 
mination of Teif and Lboi- 

4.2. Photoionization modeiiing 

Photoionization models were computed with the code nebu (e.g., 
Pequignot and Tsamis |2005) , which was also used in a com- 
prehensive study of the Sagittarius dwarf galaxy PN population 
(Dudziak et al. l2000t Pequignot et al. l2000t Zijlstra et al. l2006l ). 
Computations were done assuming that both the gas filUng fac- 
tor and the covering factor (solid angle as seen from the central 
star) were unity. The actual geometry of the nebula is not known 
due to lack of imaging and the inner radius Ri^ is a first free pa- 
rameter. A two-sector model is suggested by the coexistence of 
a very strong He ii /14686 line with [O ii] and [S ii] lines which 
implies that even though the nebula is strongly matter-bounded 
along most radial directions, some directions must be radiation- 
bounded to a large degree (One-sector models are however con- 
sidered in AppendixlA]). The premise of two-sector models com- 
prising both optically thick and thin components (advocated by, 
e.g., Clegg et al. ll987l in their analysis of the typical Galactic PN 
NGC3918) is a realistic one, given that the majority of nearby 
spatially-resolved PNe are replete with optically thick inhomo- 
geneities in the form of clumps/filaments/torii embedded in a 
more tenuous medium. The [S ii] doublet ratio can be used to 
constrain the density or, more conveniently, the thermal pressure 
Pout of the (peripheral) thick clumps. The [Ar iv] doublet ratio 
is neither suitably accurate nor sensitive enough to «e to con- 
strain the inner high-ionization region density, A^in(H) (equiva- 



0.74 ^M/Mg^1 .3 





Figure 3. The theoretical, H-burning WD tracks from VW94 are 
represented here in the log L/Lq vs. log plane. The four pan- 
els show four different metallicities, identified by the label, and 
for each metallicity tracks of different initial mass are plotted. 
A label near the knee of each track identifies its initial mass, in 
solar units. The open box represents the permitted range for the 
central star of Bennu, with parameters determined via photoion- 
ization modeling of the nebula. 

lently. Pin), which is therefore a second free parameter. Adopting 
a sufficiently flexible description for the gas distribution allows 
the exploration of a realistically large space of solutions. In prac- 
tice, the pair of (Te, «e) is obtained at each point by solving iter- 
atively the statistical balance equations, assuming a variable gas 
pressure, P, given here as a function of the radial optical depth, 
T, at 13.6 eV: 



P(t) 



Pout + Pir 



■ tan 



.log(^) 



(1) 



At the first step of the computation (t = 0) the pressure is Pin, 
while at the last step (r - oo; in practice »1) it is Pout- The re- 
quirement for a sector (Sector 2) comprising radiation-bounded 
clumps immersed in a dilute medium dictates that Pin < Pout - The 
third parameter r^- at which P(tc) = (Pin + Pout)/2 corresponds 
to the inner boundary of the generic clump. The fixed fourth pa- 
rameter, K = 30, is large enough to ensure a sharp clump bound- 
ary. For consistency, the optical depth ti of the dilute matter- 
bounded Sector 1 must be substantially less than Tq (in practice 
the condition tJti > 1.2 was adopted). 

The ionizing spectrum is described as a black body of tem- 
perature Tbb and luminosity L. According to stellar atmosphere 
models for hot stars, the 'colour temperature' of the bulk of 
the photoionizing continuum (the best equivalent black-body 
temperature) is larger than the actual effective temperature Tgff. 
When black-body spectra are used instead of genuine stellar at- 
mosphere models in computing photoionization models, a rec- 
ommended rough correction, adopted here, is: 



Figure 4. The dependence of the progenitor mass on the metal- 
licity of the best-fitting WD track is illustrated by this plot. The 
lines are least-square fits to the discrete points corresponding 
to three of the four metallicities available in VW94 tracks (see 
Fig. 3). The fits have been repeated for the two combinations of 
(T^e - AT^e; L - AL) (dashed Une) and {T^ + AT^; L + AL) (soUd 
line) allowed by our error estimates. The vertical dotted lines 
show the permitted metallicity range determined from the nebu- 
lar emission lines, and the horizontal dotted lines mark the limits 
of the permitted progenitor mass. 

Table 4. Constraints on model parameters. 



Parameter 


Constraint 


free 


standard IxlO" cm 


yVi„(H) free 


Pin < PoLtt; stdd 1200 cm-^ '[Ariv] ^4740/^47 11 -n' 


Tbb free 


tJti > 1.2; 'upper limit [0 i] /16300' 


Ti(thin) 


absolute flux /(H/3) 


Z^™^ (thick) 


[On] A3727; (/™" + ' = 1) 


^out 


ratio [Sn] ^6731/^6716 




He I ^5876 (see text) 


L 


ratio [Ar iv] /14740 / [Ar iii] /17135 


He 


He n A4686 


C 


ratio [0 ra] /14363/(/i4959+/i5007) 


N 


[N n] ^6584 


O 


[0 m] ^5007+^4959 


Ne 


[Ne m] A3869 


S 


S/Ar = 4.37 by number (solar); '[S ra] /i6312' 


Ar 


[Ar in] ^7135 + [Ar iv] ^4740 


Fe 


Fe/H = 4x10"' by number; 'upper limit [Fe vi]' 



T'eif = rBB-15kK. 



(2) 



(see, e.g., Dudziak et al. 12000b . 

The model parameters and corresponding dominant observa- 
tional and technical constraints are listed in Table |4] Some po- 
tential constraints which were dismissed on account of the large 
uncertainties they entail are given in inverted commas. Thus, we 
assumed that S/Ar was solar (4.37 by number; Lodders [20031 1. 



Table 5. Range of model parameters for 'exact' solutions 



Parameter 1-crbox 



(10" cm) 


1.0 


1.0 


1.0 


1.4 


_ 


7Vi„(H) (cm-3) 


900 


1200 


1600 


1200 


_ 


Tub (kK) 


127±5 


135±7 


145±3 


135±7 


135±10 


Teff (kK) 


112±5 


120±7 


130±3 


120±7 


120±10 


L (10'"' ergs-') 


11. ±.1 


10.±.1 


9.±.1 


11. ±.1 


10.±1. 


He/HxlO^ 


108±0 


108+0 


109±1 


108±0 


108±1 


C/HxlO^ 


197±1 


212±4 


230±4 


222±4 


214±17 


O/HxlO" 


165 ±7 


172±9 


178±4 


169±9 


170+12 


N/OxlO^ 


149±4 


147±1 


150±1 


147±1 


149±3 


Ne/OxlO'' 


142±3 


143±1 


153±2 


149±1 


145+5 


Ar/OxlO"* 


232+6 


228±5 


226±2 


228±5 


230±7 



Also, no iron lines were detected and Fe/H was given a constant 
arbitrary value. There are 1 1 (+2 implicit for S and Fe) spectro- 
scopic constraints controlling as many model parameters, leav- 
ing Tbb as a third free parameter, in addition to R^n and Mn(H) 
(see above). Gas cooling is dominated by carbon and oxygen 
line emission. For lack of measured carbon lines (no UV cov- 
erage), the carbon abundance is constrained by the energy bal- 
ance through the [O iii] ratio temperature. Given the three free 
parameters and assuming that all other constraints are system- 
atically fulfilled, for increasing L, the computed [Ar iv]/[Ar m] 
increases and He i decreases moderately, while, for increasing 
tJti, [Ar iv]/[Ar iii] increases slowly and He i decreases again, 
until a limit is reached asymptotically for tJti > 2.5-3. Thus L 
and tJt[ are primarily controlled by the ionization balances of 
argon and helium respectively. Other correspondences between 
parameters and constraints in Table |4]are straightforward. 

The domain of 'exact' solutions according to the criteria 
listed in Table |4]has been scanned. Many iterations were needed 
to fulfill all constraints simultaneously and not all combinations 
of the 3 free parameters led to a solution. Results are summa- 
rized in Table |5] Given Rin - IxlO'^cm (our 'standard' value), 
the full range of T^b leading to a solution was determined for 
the 3 given values of M„(H) (cm-^) = 900 (Pout/An = 1.35), 
1200 (standard, Pom/P\n = 1-95), and 1600 (Pout/An = 2.4). The 
ratio [Ar iv]/[Ar m] being accounted for (correct choice of L), 
a minimum Tbb is obtained when He i happens to be computed 
too large and a maximum Tbb is obtained when tJti decreases 
down to 1.2. After several trials, the resulting accessible domain 
on the (rBfi/kK, L/IO^*' erg s"') plane was found to be an elon- 
gated ellipse, whose long-axis boundaries were at (120, 11.3) 
and (150, 8.7) respectively for Mn(H) ~ 800 (the minimum value 
since Pout/An ~ 1-2), and ~ 1700 cm""*: the larger the Mn(H), the 
more compact the nebula and the smaller the corresponding L 
are. Thus, for a given Ri„, L correlates inversely with Tbb- The 
'Icr box' for the standard case, not including observational un- 
certainties, is TBE/kK = 135±10 and L/10^^*' erg s ' = 10±1. For 
Mn(H) = 1200, 'exact' models with R JIO" cm = 1.4 instead of 
1.0 were also obtained for a few Tbr's, showing that derived pa- 
rameter values were about similar to those obtained previously, 
except for a 10% upward shift of L (Col. 5, Table|5]l. Conversely, 
for Pin/lO'^cm < 1.0, the previous solutions for standard Ri„ 
essentially apply. Model outputs are displayed in Table lA.ll of 
Appendix lAl including one solution (0135, Col. 3) belonging 
to the domain of most favorable values of model parameters 
(Table |5]l and examples of models which are in some way un- 
satisfactory, including one-sector models. A more thorough dis- 
cussion of the uncertainties is also given in AppendixlBl 



Model elemental abundances are relatively stable within the 
domain of 'exact' solutions (Col. 6, Table|5]l. Model outputs (see 
Col. 3 of Table lA.ll ) suggest that any value of S/Ar between 
solar (adopted) and twice solar is possible. For higher Teff's, 
O/H can be somewhat larger and Ar/O slightly smaller (Col. 4, 
TablelAlJ- 

4.3. Parent Star 

For the central star, in the previous section we obtained best val- 
ues Teif = 120 + 10 kK and Lboi = 2600 ± 260 Lq. To determine 
an age and mass for Bennu we have compared these model es- 
timates to the H-burning white-dwarf (or PN central star) evolu- 
tionary tracks of Vassiliadis & Wood ( 119941 VW94), as shown 
by Fig. [3] As the metallicity increases, the progenitor mass of the 
best-fitting track increases. To have a quantitative estimate, for 
each metallicity the tracks were interpolated at the (log T^ff, log 
L) position of the WD, thus obtaining a mass-metallicity rela- 
tion, which is shown in Fig.|4] Due to the small number of tracks 
and the large extrapolation from the PN location to any of them, 
the Z = 0.008 point was omitted from the fit. The procedure 
was repeated for the two extremes of the (log Tg; log L) combi- 
nations allowed by the errors on the two quantities. By fitting the 
discrete points with linear regressions, one has continuous mass- 
metallicity relations, so entering with a metallicity value an es- 
timate of the progenitor mass can be obtained. If we adopt the 
widest metallicity range discussed below ([m/H] - -1.03 from 
[Ar/H] to [m/H] = -0.46 from [O/H]), then the range of pro- 
genitor mass is between 0.74 and I.3OM0. 

This low mass chimes well with the N/O ratio and He/H neb- 
ular abundance, both being sufficiently low to indicate a non- 
Type I, low central star mass status (according to e.g. the crite- 
rion of Kingsburgh & Barlow, |T994|i. From the VW94 tracks, the 
nebular age would be 23 000 yrs at the lower mass estimate, and 
7 000 yrs at the higher. To convert the progenitor mass range to 
a main-sequence (MS) age range, we looked at the isochrones 
of Pietrinferni et al. (2004) . for two values of the metallicity 
Z = 0.002 and Z = 0.008 ([m/H]= -1 and -0.4, respectively). 

The lower mass value is just below the minimum isochrone 
mass (0.77 M0), so the star might be as old as the Universe. 
To convert the higher mass value into an age, one needs to 
decide between canonical core convection, or convection with 
overshooting, which yields larger ages. In the case of canoni- 
cal convection, one obtains 2.0 or 2.7 Gyr for the two metallic- 
ities, while tracks with overshooting give 3.3 and 4.0 Gyr, re- 
spectively. Thus the overall permitted age range is greater than 
2 Gyr. 

As discussed below, it is also quite possible that the metallic- 
ity given by Ar (and S) is closer to the real PN value. In that case, 
if we take the model argon abundance, then [m/H] = -1.03+0.1, 
and by repeating the exercise above, the resulting mass range is 
0.70-1.04 Mq. Taking only isochrones at Z = 0.002, the mass 
interval is then converted into an age interval from 4.7 Gyr to the 
age of the Universe. 

5. Discussion 

5.1. Oxygen as a metallicity tracer? 

The gas-phase elemental abundance in nebulae which is the eas- 
iest to compare to stellar values is that of oxygen (i.e. [O/H]); 
[Oiii] lines are collectively strong and usually contain most 
of the flux from oxygen ions (except in cases of extremely 
metal-poor nebulae), and thus measurement and abundance- 



determination errors are small. Along with sulphur, argon and 
neon, oxygen is generally considered to be unaffected by en- 
dogenous processes within the progenitor stars. For instance, 
when comparing the oxygen abundances in PNe to those of H ii 
regions in a given galaxy the former are expected to be always 
lower (or equal) than the latter; H ii regions after all represent the 
present-day metallicity of the ISM, whereas PNe originate from 
older stellar populations when galactic metallicities are likely to 
have been smaller. But this is not always found to be the case. 

A number of studies of PNe in low-metallicity environments 
have pointed to a possible enhancement of O/H ascribed to 
'self-enrichment' processes occurring during the 3'^'^ dredge-up 
phase. PN vs. Hii region abundance comparisons reveal a pos- 
itive difference between the two at low O/H values (at about 
logO/H< -4; see e.g. Pena et al. 20071 Richer & McCall^OOT^. 
For NGC 3 109 (Pena 2007) this is about 0.3 dex, and Kniazev et 
al. ( 120071 ) adopted a similar value of 0.27 dex to 'correct' their 
PN oxygen abundance for 'self-enrichment'. Comparison be- 
tween models and observations of PNe in the Magellanic Clouds 
(Leisy & Dennefeld 120061) have also suggested that oxygen in 
PNe can be enhanced relative to the precursor metallicities, par- 
ticularly for the lower mass progenitors. These studies have in- 
dicated that argon or sulphur should be used instead of oxygen, 
as an indicator of the metallicity of the ISM from which the PNe 
formed, with argon preferred as having the more secure spectral 
determination; Pequignot et al. (120001 ) and Pequignot & Tsamis 
( I2005I I reached similar conclusions. 

It is uncertain whether Bennu is likely to contain endoge- 
nous oxygen as there are no H ii regions in Phoenix to compare 
to. Its O/H is however higher (by > 0.3 dex) than the range of 
values found in those PNe for which an oxygen enhancement 
has been suggested (although see Zijlstra et al. 120061) . In addi- 
tion, following Kniazev et al. (120071) and based on the nebula's 
Ne/O, S/O and Ar/O ratios, we should apply no correction for 
the presence of endogenous oxygen. However, the [Ar/H] abun- 
dance is lower than [O/H], so if we adopted the former as an 
indicator of the metallicity of the parent ISM, as suggested by 
Leisy & Dennefeld, then we should indeed conclude that en- 
dogenous oxygen is present in the nebula. Since with our present 
knowledge we have no exact way to decide upon this issue, in the 
rest of the discussion we will adopt a metallicity for the Phoenix 
galaxy within the range given by [O/H] and [Ar/H], and with the 
understanding that the more metal-poor end is favoured. 

5.2. The age-metallicity relation of Phoenix 

In the previous sections we concluded that the spectral prop- 
erties of the newly discovered PN in Phoenix are consistent 
with a progenitor mass of 1.0 + 0.3 and a metallicity 
[m/H] = -0.75 + 0.29 dex (or Z = 0.002 to Z = 0.008), cor- 
responding to an age of 7.9 + 5.9 Gyr Or if argon and sul- 
fur abundances are favored over oxygen, then the progenitor 
mass and age range are 0.87 ± 0.17 Mq and 9.2 ± 4.5 Gyr, for 
[m/H] = -1.03 + 0.1. 

Since this is the first direct spectroscopic determination of 
this galaxy's metallicity, it also represents the first solid con- 
straint to its age-metallicity relation (AMR). It is therefore in- 
teresting to see how this fits within our present knowledge of 
the chemical evolution of Phoenix and its star-formation history 
(SFH). 

The evidence accumulated since its discovery by Schuster & 
West (119761) and its recognition as a dwarf galaxy by Canterna & 
Flower ( |I977| l, clearly shows that Phoenix has had an extended 
star formation. The presence of an old, globular-cluster Uke pop- 



ulation in Phoenix was first established by the discovery of an 
extended horizontal branch (HE) atV ~ 23.8 (H99; M99). The 
metallicity of the intermediate-old populations was estimated 
from the color of the RGB by van de Rydt et al. ( 119911 ) and 
H99, as [Fe/H]= -2.0 dex and [Fe/H]=-1.81 ±0.10 dex, respec- 
tively. A significant spread in RGB color was also found which, 
if it was only due to a metallicity range, would correspond to a 
dispersion of about 0.5 dex. A very low metallicity at an age of 
~ 13 Gyr is then the first point in the AMR of Phoenix. 

A number of stars above the RGB tip were interpreted as 
AGB stars by H99. They are representatives of an intermediate 
age population (3 to 10 Gyr), which was estimated to be about 
30-40% of the total stellar population of the galaxy. The metal- 
licity of these stars could not be estimated, but given the RGB 
width, it could be even up to -0.3 dex, if we adopted an upper 
Umit of 3o- (i.e. -1 .8 + 3 x 0.5 dex). 

The presence of a young stellar population was established 
quite early (Canterna & Flower [T977l Ortolani & Gratton [T988] ), 
and together with the possible association of an H i cloud, lead 
to a morphological classification for the Phoenix as intermedi- 
ate between dSph and dirr galaxies (e.g.. Young et al. 2007; 
Carignan, Demers, & Cote, 119911 ). According to H99 the star 
formation episode started at least 0.6 Gyr ago, and lasted until 
10^ yr ago. Since the galaxy has no Hii regions, and since its 
main-sequence stars are too faint, no spectroscopic determina- 
tion exists for the current metallicity of Phoenix. 

The current metallicity and future evolution of Phoenix de- 
pend on the presence of gas associated with the galaxy. The 
recent measurement of the optical radial velocity of Phoenix 
(Irwin & Tolstoy 120021 ) has shown that it is the same as that of 
'cloud A' of Young and Lo ( 119971 ). confirming a long-suspected 
association (St-Germain et al. [1999). This means that the Hi 
mass is ~ lO^M©, and as shown in H99, it can be explained 
by mass lost from RGB and AGB stars, and PNe, over the last 
~ 2 Gyr And according to Young et al. (2007) the gas will 
not be able to escape the galaxy, so SF might start again in the 
future. However currently SF from molecular gas is inhibited, 
since no CO was detected by Buyle et al. (120061 ). and Jackson 
et al. (120061 ) found no diffuse 8 yum emission from dust. The 
(unknown) metallicity of Phoenix youngest stars therefore rep- 
resents the present end point of the galaxy's chemical evolution. 

It is also important to note that age gradients exist throughout 
the galaxy, with all young stars concentrated in an inner com- 
ponent, while the old population is more extended and oriented 
N-S (M99; H99). The fraction of intermediate-age stars over old 
stars also decreases going from the central to the outer regions 
(Hidalgo et al. 2007). Moreover, additional support for an ex- 
tended SFH is given by the study of the variable-star content of 
Gallart et al. ( 120041 ). They find both anomalous Cepheids and 
short-period classical Cepheids, which can be explained if the 
metallicity has been low ([Fe/H]= - 1 .3) for most of the galaxy's 
lifetime. 

The star formation history for the central part of Phoenix 
has been derived by Holtzman et al. (120001 HOO) and Young et 
al. (2007), using a WFPC2 CMD reaching the oldest main se- 
quence turnoffs. Hidalgo et al. (2007) present a determination of 
both the star formation rate as a function of time (SFR(t)) and 
the metallicity as a function of time (Z(t)) for the same field and 
an outer one. All these authors agree on that Phoenix has had 
an almost continuous SFH, and with a roughly constant star for- 
mation rate, somewhat decreasing toward the present time. The 
star formation rate at intermediate ages (^ 6-2 Gyr) seems to 
have been somewhat larger than immediately before and after, 
specially in the central part of the galaxy where Bennu resides. 



Integrating their age distribution one gets that 15% of the stars 
have ages between 1.5 and 3 Gyr, and another 15% have ages 
between 3 and 10 Gyr, which is consistent with the fraction of 
intermediate-age stars estimated by H99 (see above). Ahhough 
HOO do not give a real AMR, from their 'population box' one can 
see that the [Fe/H] increases from ~ -2 at the earliest stage up to 

-(1.5-1) at ~ 3 Gyr, then reaches 1-0 at ~ 1 Gyr and stays 

constant thereafter. In the case of H07, they predict quite low 
metallicities up to the current epoch. The maximum is reached 
ca. 1 Gyr ago, and it is less than -1 in [m/H]. 

Turning now back to Bennu, we can place its representa- 
tive point in the AMR relations of HOO and H07, which are 
the only two studies making quantitative derivations. In the age 
range of Bennu's progenitor, most of the simulated stars of 

HOO have metallicities around [m/H] 1, but there is also 

a smaller group of objects with metallicities centered around 
[m/H] - -0.5. Adopting instead the age range suggested by ar- 
gon and sulfur abundances, the metallicities of HOO are around 

[m/H] 1 .5 with a tail up to [m/H] = - 1 . Therefore the metal- 

licity we measure tends to be higher than that of HOO in the same 
age range, but still compatible within the uncertainties. A metal- 
licity lower than ours is also predicted by the simulations of H07, 
whose AMR generally stays more metal-poor than that of HOO. 
As it was recalled above, even the peak of the relation happens 
at a metallicity lower than that of the PN. While an independent 
study would be needed to resolve this issue, we can conclude 
that existing AMRs predict, for the age of Bennu, a metallicity 
that is lower than the one we measure, particularly so for the H07 
study. 

5.3. Phoenix in tlie context of Local Group galaxies 
5.3.1. The number-luminosity relation 

A recent review of the PN population in Local Group (LG) 
galaxies can be found in Corradi & Magrini (120061 CM06), 
where they list 20 galaxies containing at least 1 nebula. These 
objects belong mostly to the northern sky, where the most sys- 
tematic searches are on-going. There is a very good linear cor- 
relation between the V luminosity of a galaxy and the number 
of PN candidates, which predicts that below Ly = 10^ Lq the 
probability of finding a PN is very low. This means that in prin- 
ciple Phoenix, which has aLy - 10^ Lq, should have no PNe. 
However the number of nebulae also depends on the SFH of a 
galaxy, and in particular an enhancement at intermediate ages 
can also increase the number of detected PNe (see Buzzoni et 
al. 120061 for a thorough treatment of the luminosity-specific fre- 
quency of PNe). Some evidence of this possibility is given in 
CM06, where it is shown that galaxies with enhanced SFH in 
the age interval 2 to 8 Gyr have indeed a larger number of PNe 
(for example NGC 205 and the SMC). The presence of a PN in 
Phoenix despite its low luminosity could therefore be a sign of 
substantial SFH at intermediate ages. This is nicely consistent 
with the large number of AGB stars detected by H99, who, as 
recalled above, estimated that ~ 37 ± 12% of the total popula- 
tion should consist of stars with ages of 3 to 10 Gyr. The en- 
hanced SFR between 2 and 3 Gyr found by HOO (for a central 
field very close to Bennu) is also in line with this finding. An 
additional possibility, discussed in the next section, is that the 
luminosity of Phoenix is underestimated, either because of mea- 
surement uncertainties, or because of mass being stripped from 
the galaxy by the tidal field of the Milky Way. A combination of 
these two effects might be able to produce a luminosity which is 
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Figure 5. The luminosity-metallicity relation for Local Group 
galaxies, with data from Corradi (priv. comm.); the [Fe/H] val- 
ues come from stellar results (old-intermediate populations) 
and [O/H] is from nebular abundances (PNe and Hii regions; 
intermediate-young populations). My were taken from Mateo 
(119981 ). The dashed line represents an unweighted linear fit, 
while the dotted lines represent the +lcr dispersion around the 
average fine. 



a factor of 10 larger than what is currently published, and so to 
reconcile Phoenix with the number-luminosity relation. 

5.3.2. The luminosity-metallicity relation 

In Fig.|5]we plot the iron and oxygen abundances vs. the absolute 
V luminosity for a sample of Local Group galaxies (R. Corradi, 
priv. comm.). The asterisk represents Phoenix, with [Fe/H] from 
H99 and 12 + log (O/H) from this paper. Its error bars are the 1- 
sigma dispersion for the iron abundance given by the spread in 
color of the RGB, while for the oxygen abundance the bar rep- 
resents our permitted range. The luminosity has been computed 
by taking the apparent luminosity V - 13.2 from de Vaucouleurs 
& Longo (119881) and converting it into My - -10 ±0.13 
with the distance modulus of H99 based of the HB luminosity: 
(m - M)o - 23.21 ± 0.08. An error of 0.1 mag was assumed for 
the catalog luminosity. Since the relation is steeper for [Fe/H], 
the figure might be telling us that lower-mass galaxies undergo a 
comparatively larger metal enrichment than higher-mass galax- 
ies, when going from intermediate-old populations to younger 
populations (note that the range of the y-axis is the same in both 
panels). In a simple closed-box picture of chemical evolution, 
this would mean that smaller galaxies have been able to convert 
larger fractions of their gas into stars. However the figure might 
also be telling us that [O/Fe] is enhanced in low-mass galaxies, 
as discussed in Sect. 15. II The scatter of the relations is probably 
due to a combination of measurement errors, intrinsic abundance 
scatter, and age differences of the metallicity tracers. Indeed old 
to intermediate stars contribute to the [Fe/H] values estimated 



with the color of the RGB, and intermediate to young popula- 
tions contribute to the [O/H] values from PNe and H n regions. 

Figure |5] shows that the value of [O/H] of Phoenix is more 
than 1-sigma above the average, for its luminosity. There are 
several possible explanations for this fact. First, it is possible 
that the PN's oxygen abundance would lead us to overestimate, 
as discussed in Sect. 15.11 the metallicity of the galaxy. Indeed 
recent theoretical AGB models by Marigo ( 1200 It MOl) predict 
positive oxygen yields for 0.8 < M/Mq < 3.5, with a strong de- 
pendence on metallicity. This is due to the fact that, as metallicity 
decreases, the thermal-pulse AGB phase increases its duration. 
Therefore it allows more dredge-up episodes, which also happen 
to be more efficient. As we found above, the progenitor star of 
Bennu is in the range 1 < M/Mq < 1.6 and so it is formally 
in the range where oxygen enhancement is expected. However 
in MOl models the yield has a maximum for masses between 
2Mq and SM© and it becomes very modest below ~ 1.3Mq. 
On the other hand the minimum metallicity in MOl models is 
[m/H] - -0.7, and we might expect a larger oxygen enhance- 
ment if the metallicity of Bennu is closer to -1 as suggested by 
the argon abundance. Lacking a custom theoretical model con- 
structed for the Phoenix PN progenitor, the conclusion is that an 
endogenous production of oxygen might explain the position of 
Phoenix in the luminosity-metallicity relation (see also Magrini 
et al. 12005 I for an extensive discussion of oxygen enhancement in 
PNe). Note also that, taking argon as a metallicity indicator (i.e. 
the lowest point of the error bar in the figure), the abundance 
would be within the r.m.s. dispersion of the relation. 

It is also possible that the progenitor was born within an ISM 
recently enriched by Type II SNe. But perhaps the problem lies 
in the luminosity of the galaxy. The V luminosity is flagged as 
'Low quality data' in the NASA/IPAC Extragalactic Database, 
and for example the R luminosity from Lauberts and Valentijn 
( II989I I is 1 magnitude brighter. A luminosity of ~ - 11 would 
bring the galaxy within the 1-sigma dispersion of the relation. 
And finally another possibility is that the current mass of the 
galaxy is not representative of the mass at the time when the PN 
was born. In fact M99 find that the radial profile of Phoenix can 
be well fit by a King model, which suggests a tidal truncation 
by the gravitational field of the Milky Way. It also means that 
part of the galaxy mass was lost due to dynamical evolution, and 
hence that the luminosity of Phoenix was larger in the past. 

6. Summary and conclusions 

A planetary nebula was recently discovered in the Phoenix dwarf 
galaxy, the first ever found in this stellar system. In this paper 
we presented our follow-up spectroscopic data obtained at the 
ESOA^LT with FORSl. With a total exposure time of 1.4 hours 
we measured emission line fluxes down to 9x 10"'^ erg cm"^ s"' . 
Such high sensitivity allowed to detect, together with the Balmer 
series, all lines of oxygen, neon, sulfur, helium, argon, and nitro- 
gen above the quoted threshold. To calculate the abundances of 
these elements, both the empirical method, and photoionization 
modeling were employed. The element-to-hydrogen abundance 
ratios were found to be consistent with each other for the two 
methods, but the second method yields more reasonable lumi- 
nosity and temperature for the central star. The oxygen abun- 
dance is larger than that of argon and sulfur. This might be due to 
endogenous production of oxygen in the progenitor, or it might 
be caused by a truly higher oxygen abundance of the ISM where 
the progenitor was born. To decide between the two hypotheses 
an independent measurement of the abundance would be needed, 
which however does not exist. Therefore we based our discus- 



sion on a metallicity range comprised between [m/H]= -1.03 
(sulfur and argon abundances) and [m/H]= -0.46 (oxygen abun- 
dance). Our conclusions would be essentially the same if we re- 
stricted the metallicity range to that allowed by argon and sulfur 
alone. 

Using VW04 tracks and Pietrinferni et al. ( 120041 ) isochrones, 
we found that the progenitor star of Bennu should have an age 
comprised between 2.0 and 13.7 Gyr. The more restrictive low 
metallicity would give an age 9.2 + 4.5 Gyr So despite a con- 
siderable uncertainty, the Phoenix PN allows to put a constrain 
on the galaxy's age-metallicity relation. This shows that exist- 
ing AMRs underestimate the metallicity at intermediate ages, by 
as much as ~ 0.6 dex, even when adopting the most restricted 
[m/H] range. 

Within the general picture of LG dwarf galaxies, the pres- 
ence of a PN in Phoenix is unexpected given the galaxy's low 
luminosity (Corradi & Magrini r2006b . Moreover, an extrapola- 
tion of a linear L-Z relation to the luminosity of Phoenix would 
suggest that the metallicity of the galaxy is larger than that of 
other LG galaxies of comparable luminosity (at an intermedi- 
ate age). On the other hand there are no other galaxies at such 
low luminosity with measured nebular oxygen abundance, so it 
might be possible that the relation simply deviates from linear- 
ity for the lowest-mass galaxies. Another possibility is that the 
problem resides in the luminosity. The galaxy might be more lu- 
minous than quoted in the literature, so a modern measurement 
of the integrated light would be very desirable. It is also very 
likely that Phoenix lost a fraction of its mass through tidal inter- 
action with the gravitational potential of the Milky Way (M99), 
so its representative point in the L-Z relation should me moved to 
higher values of L. This would also help with the PN-number vs. 
luminosity relation, thereby explaining the presence of Bennu. 

In conclusion, this work demonstrates that the detection of 
even a single PN in a dwarf galaxy, can provide essential infor- 
mation on its chemical and dynamical evolution, and it can lead 
to a better understanding of its past star formation history and 
mass build-up. 
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Appendix A: Photoionization model outputs 

Two-sector models ('D') using standard assumptions (Table |4]l 
and one-sector models ('S'), labelled by T^b in kK, are pre- 
sented in Table lA.fl All but one (D135) fall outside the domain 
where 'exact' solutions can be found (Sect. 14.21 Table |5]l, and 
'best fit' values are adopted for He i, [Ar m], [Ar iv] and even 
He II. D105 cannot account for He ii for any He/H and ti, and the 
computed He i is unacceptably large. D160 illustrates the max- 
imum Tbb obtained when the computed He i /15876 intensity is 
allowed to depart from observation by at most 33%. For even 
larger Tbb, He i is underestimated and the ratio [Ar iv]/[Ar iii] is 
not accounted for either Also, the optically thick sector can de- 
velop a warm neutral zone beyond the ionization front that emits 
(unobserved) [O i] ^6300. 

Computed [S ii] and [S iii] intensities for 0135 suggest that 
S/Ar may be up to twice larger than the adopted solar value, but 
this is not judged significant, given the weakness of observed 
[S III] i6312. From the strongest predicted iron line, [Fe vi] 
/15145, an upper limit to Fe/H is twice the adopted value, so that 
Fe/Ar is less than 0.25 solar; as in usual PNe, most of the iron is 
likely to be locked into dust grains. In Table lA. fl the ratio [Ar iv] 
/1471 1H-//14740 increases with Tbe due to the growing contribu- 
tion of [Ne iv]. 



For single-sector uniformly matter-bounded models to si- 
multaneously account for He ii and [O ii] requires, at the same 
time, a high Tbb, a low L and small ionization parameter (large 
Rin and/or large A^in). These models typically fail to reproduce 
the argon ionization balance, underestimating the [Ar iv]/[Ar iii] 
ratio by factors of >3, and imply large He/H (in the models dis- 
played, He/H is too small). These models are rejected, as well 
as their low value of O/H. Any positive detection of [O i] would 
further invalidate the single-sector assumption. 

Appendix B: The central star parameters 

Uncertainties affecting the model geometry and the spectro- 
scopic constraints can enlarge the domain of acceptable solu- 
tions in the (Tbb, L) plane. Uncertainties on L, related to Rm and 
to [Ar iv]/[Ar iii], amount to a factor of 1.6 and are therefore 
not too influential in the interpretation of the stellar parameters 
(Sect. |5]l. This factor, obtained assuming a factor 1.45 uncer- 
tainty on [Ar iv]/[Ar iii], allows for both observational and the- 
oretical uncertainties (the recombination coefficients for argon 
are inaccurate, e.g., Dudziak et al. 120001) . The range of accessi- 
ble TsB is controlled by the assumed value of A^in(H) (Table [S]) 
and, particularly, the intensity of He i /15876. The computed He i 
increases rapidly towards low values of T^b (alias Tetj) because 
He/H must increase to account for He ii (see model D105, Col. 2 
of Table lA. ll . In addition, the assumption of a strict black body 
for the central star, useful in accounting for the strong He ii 
line, is more likely to break down for lower Teff's due to the 
occurrence of a discontinuity at the ionization limit of He^ (e.g., 
Rauch l2003l l. thus exacerbating the problem. By itself, the in- 
crease of He/H strongly militates against too small values of 
Tbb- Thus, a lower limit to Tbb is likely close to 120kK, cor- 
responding to a minimum Tgff ~ 105 kK (Eq. 2). On the other 
hand, trial computations showed that dividing the present He i 
line intensity by 1.33, models could be obtained up to Tbb = 
160 kK (D160, Col. 4 of Table IaTTT i. equivalent to a maximum 
Tgff ~ 145 kK. In fact, a few models re-converged using stellar 
atmosphere models (Rauch 120031 see Fig. 22 in Pequignot & 
Tsamis l2005b showing that a more nearly correct maximum Teff 
is almost 150 kK. More accurate He i line intensities are needed 
before the range of possible Teff's can be narrowed. Note that 
the derived standard He/H {- 0.108) is relatively large, suggest- 
ing that either the observed He ii line intensity is slightly over- 
estimated or some basic model assumption is lacking. To our 
knowledge, only in a situation of chemical inhomogeneity of the 
PN could the derived He/H be lowered (Pequignot et al. 2002). 
Indications in favour of such a possibility are: (i) Hy is ob- 
served to be stronger than expected, and (ii) the [S n] +C m /14069 
and C iv+[Fe iii] /14658 blends are severely underestimated in 
all models (Table lA.lb . These spectral regions include C and 
O recombination multiplets which in many PNe can only be 
accounted for with models that harbour H-deficient inclusions 
(Pequignot et al. 2002). The current observations of Bennu are 
not sufliciently deep to allow us to pursue this possibiUty further. 



Table A.l. Model parameters and results. 



Model" 




D105 


D135 


D160 


S180 


S310 






2423 


2550 


2420 


2080 


2084 


Ri„ (10" cm) 




1.00 


1.00 


1.00 


2.58 


1.82 


RUthin) (10" 


cm) 


2.08 


2.15 


2.14 


2.97 


2.06 


/?out(thick)(10"cm) 


2.91 


2.98 


2.80 


_ 


_ 


(thick) 




0.027 


0.013 


0.012 


_ 


_ 


A'i„(H) (cm-3) 




1200 


1200 


1200 


1000 


1600 


<;7e> (cm^^^) 




1942 


1860 


1870 


1590 


2940 


<r(H+)> (10'' K) 


1.752 


1.780 


1.781 


1.730 


1.669 






1.00 


1.75 


1.65 


0.90 


0.90 


T(thin) (13.6eV) 


0.50 


0.79 


1.09 


1.37 


3.40 


HyS(thin) fraction 


0.73 


0.87 


0.89 


1.00 


1.00 


M,o„(10-'Mo) 




6.83 


7.02 


6.85 


6.11 


3.37 


He/H 




0.122 


0.108 


0.105 


0.120 


0.110 


C/H(xlO-^) 




13.0 


22.0 


26.5 


11.7 


25.0 


N/H (xlO-^) 




1.30 


2.56 


3.10 


1.42 


1.24 


O/H(xl0-5) 




8.35 


17.5 


22.0 


7.72 


8.26 


Ne/H(xlO-') 




1.22 


2.57 


3.49 


1.29 


1.50 


S/H(xlO-'') 




14.9 


17.3 


19.7 


12.7 


11.4 


Ar/H(xlO"') 




3.40 


3.96 


4.55 


2.90 


2.60 


Fe/H(xlO-'') 




4.00 


4.00 


4.00 


4.00 


4.00 


Lines 


Jobs 


Predicted intensity / Observed intensity 


Ha 6563 


284. 


0.99 


0.99 


1.00 


1.04 


1.06 


Hy 4340 


50.4 


0.94 


0.94 


0.94 


0.94 


0.94 


H(5 4101 


27.7 


0.96 


0.97 


0.97 


1.00 


1.01 


He I 5876 


2.8 


2.56 


1.00 


0.75 


1.13 


0.74 


He II 4686 


98.1 


0.88 


1.00 


1.00 


0.98 


0.91 


C III] 1908 


— 


(591) 


(573) 


(582) 


(625) 


(1267) 


C IV 1549 




(1239) 


(1864) 


(1978) 


(667) 


(835) 


C IV 4658'- 


3.0: 


0.12 


0.18 


0.22 


0.12 


0.09 


[N n] 6584 


9.4 


1.00 


1.00 


1.00 


1.00 


1.00 


[O I] 6300 


<2. 


0.64 


0.74 


0.77 


0.03 


0.03 


[O II] 3727 


26.8 


1.00 


1.00 


1.00 


1.00 


1.00 


[O III] 4363 


8.3 


1.00 


1.00 


1.00 


1.00 


1.00 


[O ni] 5007+ 


451. 


1.00 


1.00 


1.00 


1.00 


1.00 


[O iv] 26^m 




(229) 


(551) 


(619) 


(151) 


(107) 


[Ne III] 3868 


28.0 


1.00 


1.00 


1.00 


1.00 


1.00 


[S ii] 4069" 


3.7: 


0.31 


0.22 


0.24 


0.25 


0.21 


[S ii] 6716 


1.7: 


0.60 


0.40 


0.39 


0.50 


0.41 


[S ii] 6731 


2.6: 


0.63 


0.41 


0.40 


0.49 


0.41 


[S m] 6312 


1.1: 


0.94 


0.71 


0.74 


1.29 


1.11 


[Arm] 7136 


2.7 


1.38 


1.00 


1.00 


1.63 


1.45 


[Ariv] 47ir 


5.4 


0.76 


0.95 


1.02 


0.41 


0.29 


[Ar iv] 4740 


2.7 


0.93 


1.00 


1.00 


0.47 


0.34 


[Fe vi] 5145 


<3. 


0.33 


0.48 


0.42 


0.26 


0.20 



[°] D and S refer to 2- and 1-sector models respectively. Models are 
labelled by Tbb in kK. 

[*] '<' upper limit and ':' large uncertainty. Unobserved lines: pre- 
dicted intensities are given in parentheses in units of H/? = 100. 

['■] Blends: C iv 4658 + [Fe m] 4658 + C m 4649 + O n 4651; [S n] 
4069 + [Fev] 4071 +Cin4069; [Ariv] 4711 + He 1 4713 + [Neiv] 
4714+16. 



